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Summary 

The  effective  use  of  mobile  sensors  (UAVs  for  example,  but  not  exclusively)  to 
investigate  own  and  opponent  force  status  in  geographical  regions  requires  planning  to 
compensate  for  their  limitations.  These  limitations  include  finite  endurance,  and  realistic 
mission  unreliability:  the  failure  propensities  of  platform  and  sensing  packages,  and  the 
consequent  need  for  maintenance  and  logistics  support. 

This  paper  supplies  analytical  (1 mathematical )  models  (AMM)  to  assist  planners,  and 
the  acquirers  and  operational  testers  of  mobile  sensor  assets.  The  results  are  formulas 
(that  can  be  quickly  evaluated  numerically)  for  the  expected  time  on  station  (fraction  of 
time  a  region  has  sensor  presence)  of  several  cooperatively  operating  systems  supported 
by  several  maintenance  facilities.  These  results  can  be  used  to  estimate  the  number  of 
sensors  of  a  particular  type  needed  to  cover  a  region  adequately. 

Limited  evidence  indicates  usefully  satisfactory  agreement  between  analytical  model 
estimates  and  those  from  more  detailed  and  elaborate  Monte  Carlo  simulations  (MCS). 


Even  with  modem  fast  computers  the  exercise  of  AMMs  is  more  rapid  than  of  MCSs,  so 
initial  exploratory  investigation  by  AMMs  is  indicated. 


1.  A  Generic  Situation 

A  platform  carrying  sensing  equipment,  e.g.  JSTARS,  AWACS,  or  an  RPV/UAV  is 
used  to  carry  out  surveillance  of  a  region.  The  platform  has  finite  endurance:  it  requires  a 
time  T  to  transit  outbound  from  home  base  (airbase  on  land.  Navy  carrier)  to  the  target 
region,  then  can  remain  on-station  over  that  region,  performing  surveillance,  for  time  S, 
and  finally  transit  back  again  for  time  T.  Once  back  on  base  it  spends  a  time  D  for 
refueling,  maintenance,  etc. 

Often  several  such  platforms  must  operate  cooperatively  in  order  to  achieve  adequate 
coverage  of  a  relatively  distant  region.  These  may  be  serviced  at  home  base  by  several 
maintenance  facilities;  in  some  cases  there  are  fewer  maintenance  “paths”  than  platforms, 
in  which  case  queues  can  form  and  the  coverage  correspondingly  degenerates.  We 
provide  models  that  describe  and  predict  the  long-run  consequences  in  such  cases. 

1.1  A  Basic  Formula:  Consider  that  a  cycle  begins  each  time  the  UAV  takes  off  from  its 
home  base,  and  ends  when  it  is  again  ready  to  do  so.  Let  Y(t)  =  1  if  the  platform  is  on 
station  and  Y(t )  =  0  otherwise.  { Y(t);  t  >  0}  is  a  degenerate  alternating  renewal  process. 
The  long-run  proportion  of  time  the  platform  is  on  station  if  no  systems  fail  is  seen  to  be 

lim  -  f  Y(i)dt  =  n  = - - - . 

s->cos{  2  T  +  S  +  D 


(1.1) 


0 


However,  it  is  realistic,  and  essential,  to  account  for  failures  during  the  mission,  and  their 
effect  on  time  on  station.  Failures,  and  their  repair  times,  make  the  on-station  availability 
predicted  by  (1.1)  degenerate,  and  our  analytical  models  efficiently  show  how  this 
occurs,  although  approximately.  The  Monte  Carlo  simulations  of  Pendergast  and 
Stoneman  (1998)  are  potentially  more  versatile  and  accurate,  but  require  more  setup  and 
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computer  time.  The  analytical  models  (here)  are  convenient  and  adequate  for 
explorations. 

2.  Single  System  Models 

VARIATION  1:  Assume  mission-affecting  failures  (MAF)  occur  according  to  a  Poisson 
process  with  rate  A.  Once  a  single  MAF  occurs,  the  UAV  immediately  returns  to  base 
(crashes  are  not  modeled  here).  If  additional  failures  occur  when  returning  to  base,  these 
influence  the  total  repair  time.  Model  the  repair  times  as  independent  and  random,  with 
mean  \/fi  and  let  C  be  the  duration  of  a  cycle:  total  time  from  one  takeoff  to  the  next. 

Let  X  be  the  time  until  failure.  Let  A  be  a  random  on-station  time,  and  let  B  be  the 
time  from  platform  take-off  until  return  to  base,  then  it  is  shown  in  Appendix  I  that 


(2.1) 

VARIATION  2:  Suppose  the  need  for  unscheduled  maintenance  activities  (UMA)  occur 
according  to  a  Poisson  process  with  rate  Ou  during  the  platform  flight.  These  (UMA)  do 
not  shorten  the  flight  but  each  requires  a  random  time  to  perform  with  mean  1//3 u  hours. 
In  addition  there  is  one  Administrative  logistic  down  time  (ALDT)  with  mean  1/Pa. 
Assume  that  scheduled  maintenance  activities  (SMA)  occur  at  a  rate  of  Os  per  flight  hour. 
Each  scheduled  activity  requires  l/Ps  hours.  Model  E[D]  as  follows 

E[D]=  a,jE{B}  -U  g5E[E]  j-+j- 

expected  ^  approx. 

number  of  number  of 

UMA'sto  SMA' s  to 

occur  in  a  cycle  occur  in  a  cycle 
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Numerical  Example. 

Assume  MTBMAF=25  hours.  There  is  scheduled  maintenance  of  7  hours  at  50  hr. 
intervals.  UMAs  occur  at  a  rate  of  A  =  1/5  per  hour  and  each  requires  l//3y  =1.9  hours  of 
maintenance;  the  mean  ALDT  time  1/(5a  =  0.5.  The  mean  repair  time  for  a  MAF  is  1.9 
hours  +0.5  ALDT.  The  total  time  the  platform  can  be  in  the  air  is  20  hours.  T  is  the 
ingress/egress  time. 


Fraction  of  Time  at  Least  One  UAV  is  On  Station 
_ ETOS  [Estimated  Time  on  Station] _ 


mm 

Analytic  ETOS 

Simulation  ETOS 
(95%  Confidence  Interval) 

1 

0.52 

0.5200-0.5332 

2 

0.45 

0.4389-0.4511 

3 

0.38 

0.3697  -  0.3849 

4 

0.32 

0.3121-0.3284 

5 

0.26 

0.2526  -  0.2667 

6 

0.21 

0.1978-0.2094 

3.  Many-System  Models:  Ample  Service 

VARIATION  3:  Assume  there  are  N  platforms.  Assume  there  is  ample  service  capacity,  here 
N  <  s,  the  number  of  maintenance  paths  or  servers  so  that  no  platform  needs  to  wait  for 
maintenance  or  repair.  The  platforms  are  launched  independently  of  one  another  when 
they  are  available;  a  pessimistic  assumption.  Then  an  approximation  for  the  long-run 
proportion  of  time  the  region  is  under  surveillance  by  at  least  one  platform  is  1  -  (1  -  n)N. 
For  the  parameters  of  the  previous  numerical  example,  wherein  N  =  4,  the  long-run 
proportion  of  time  the  region  is  under  surveillance  by  at  least  one  platform  is 
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Fraction  of  Time  at  Least  One  UAV  is  On  Station 


T 

Analytical  Mode! 
Fraction  of  Time 

5>5^  - 

i  Cbnfideiice  IntervaJl 

1 

0.95 

0.9565-0.9642 

2 

0.91 

0.9161-0.9247  1 

3 

r  085 

0.8694-0.8946 

4 

0.79 

0.8279-0.8451 

5 

0.70 

0.7535-0.7949 

6 

0.60 

0.6921-0.7093 

The  simulation  results  are  from  the  MASS  model;  Pendergast  and  Stoneman  (1998). 

The  analytical  model  gives  smaller  fractions  of  time  the  region  is  covered.  This 
behavior  is  because  the  simulation  includes  scheduling  of  the  UAVs,  which  the  analytical 
model  does  not  include. 


4.  A  Model  With  Limited  Service:  Markovian  Approximation 

Suppose  there  are  s  maintenance  paths:  servers  in  queuing  theory  parlance.  Let  the 
times  between  maintenance/repair  for  a  UAV  be  iid  exponential  with  mean  1  /y,  and 
assume  the  times  for  maintenance/repair  are  iid  exponential  with  1/p.  Let  X(t)  be  the 
number  of  UAVs  down  (receiving  or  awaiting  maitenance/repair)  at  time  t. 


Model:  (X(f);  t  >  0}  is  a  continuous  time  Markov  chain  with  limiting  distribution 


V(x)  =  lim P{ X(t)  =  x}  =  fl- 


i=o  min(i  +  l,s)p 


¥(0) 


(4.1) 


for  x  —  1,.  * N  where  N  is  the  number  of  UAVs.  y/(0)  is  found  by  solving 

N 

1  =  2>(0  (4.2) 

t  t  _  £[A]  .  . 

^  K v  ~  £[fi]  1)6  30  aPProximation  for  the  fraction  of  time  a  flying  UAV  is  on 


station.  Model  the  long-run  fraction  of  time  the  region  is  covered  by  a  UAV  by 
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(4.3) 


~  0  Kv) 

r=0 


Numerical  Example 

Take  l/y=£[5]  and  1/p  =  £[7)]+(y/y+— -  .  For  the  parameter  values  of 

v  ^  2/  *4“  u  / 

the  previous  numerical  example  with  4  UAVs,  7  =  3  and  a  variable  number  of  servers: 


Fraction  of  Time  at  Least  One  UAV  is  On  Station 


Number 

of 

Servers 

Model 

Fraction  of  Time 
Area  Covered 

Simulation 

95%  Confidence 
Interval 

1 

0.62 

0.8299-0.8601 

2 

0.80 

0.8754-0.8947 

3 

0.85 

10.8713- 0.8951 

4 

0.85 

0.8750-0.8891 

If  there  are  4  UAVs  ,  one  maintenance  path  and  variable  ingress/egress  time 


Fraction  of  Time  at  Least  One  UAV  is  On  Station 


7 

Model 

Fraction  of  Time 
Area  Covered 

:  >  Simulation 
95%  Confidence 
Interval 

1 

0.73 

0.9529  -  0.9604 

2 

0.68 

0.8977  -  0.9194 

3 

0.62 

0.8299-0.8601 

D 

0.55 

0.7678  -  0.7916 

5 

0.48 

0.6540  -  0.6836 

6 

0.40 

0.5230  -  0.5493 

5.  A  Model  with  Scheduling:  Replacement  UAVs  Always  Available 

In  this  section  we  describe  a  model  that  accounts  for  scheduling  the  UAVs  so  that  one 
is  always  on  station.  Initially  assume  that  there  are  always  UAVs  to  replace  on-station 
UAVs.  Let  7  be  the  ingress/egress  time  and  S  be  the  time  on  station.  We  assume  27  <  S. 

Assume  a  UAV  has  just  arrived  on  station  at  time  0.  Its  replacement  must  take  off  at 
time  5  -  7  in  order  to  relieve  it  when  the  first  UAV  starts  its  trip  back  to  base.  The  egress 
time  is  7. 
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Assume  a  UAV  experiences  a  mission  affecting  failure  after  an  exponential  length  of 
time  having  mean  1/A. 

Assume  a  replacement  UAV  is  launched  T  time  units  before  the  on-station  UAV  is  to 
start  to  return  to  base  or  upon  failure  of  the  on-station  UAV  whichever  comes  first.  Let  L 
be  the  length  of  time  a  replacement  UAV  is  late  to  the  region;  during  this  time  the  region 
is  not  covered. 


£[L]  = 


replacement 
docs  not  fail 
on  ingress 


replacement 
fails 


x  +  E[L] 


new 

placement 


dx 


(5.1) 


Thus, 

^[i-e-xr -(xry"]  1 

£[£]  =  — - JV - i[^-l -xr\  (5.2) 

Notice  that  the  expected  lateness  increases  exponentially  with  the  transit  time  to  station, 
which  argues  for  getting  closer  to  the  region. 

Consider  a  cycle  to  start  when  a  new  UAV  arrives  on  station.  The  expected  length  of 
a  cycle  is 

E[C]  =  e~^  [S  +  E[L]] + [x  +  E[L]  +  T]+  fje’^dx  [S  +  £[Lj  (5.3) 

on-station  on-station  UAV  v  <✓  ■' 1  * 

UAV  docs  not  fails  before  on-station  UAV  fails 

fail  on  station  replacement  after  replacement 

is  launched  is  launched 


=  E[L]  +  e~x{s~T)S  +  -fl - e~x{s~T)  - (A(S - T))e~x(S~T] ]  +  r[l- e~x[s'T) ]  (5.4) 

A 

=  £[£,]+-[ l-e~x{5~T)  -  A(S - T)e~X{S~T)]  +  e~X(S'T]S  +  r[l - e~x(s~T)]  (5.5) 

A 
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Let  V  be  time  the  region  is  not  covered  during  a  cycle. 


E[V]  =  (T+ E[L ]) (l - e~x(S-T))  +  e~x(S-T)fo  Xe^ [T-x  +  E[L]]dx 

^  V  J  '  - 11  1  v  * 

probability  on-sUtion  on-station  UAV  fails 
UAV  fails  before  after  replacement  takes  off 
replacement  takes  off 


-XT 


=  (T  +  E[L])(l-e'HS-T)) 

+ c--‘<5-r)|[r  +  £[L]][l  -e-*]-  j[l  -  „ 
The  long  run  fraction  of  time  the  region  is  not  covered  is 

E[C] 


(5.6) 


(5.7) 


6.  A  Model  with  Scheduling:  Limited  Repairs 

Let  Xi  be  the  number  of  failures  during  a  cycle  for  the  ingressing  UAV.  If  there  is 
always  a  replacement  UAV  available,  then 


£[X,]  =  (l-c-”')+(l-^-”')£[X,] 

(6.1) 

ElX,]  =  ^P-  =  eXT -1. 

(6.2) 

Let  Xo  be  the  number  of  failures  for  the  UAV  on  station. 

E[X0]  =  (l -  e-A(5+T)) 

(6.3) 

where  T  is  the  egress  time. 

The  long  run  average  failure  rate  is 

,  E[X0]+E[X,] 

°~~~  £[C]  • 

(6.4) 

Consider  a  queue  with  arrival  rate  Ao  and  a  maximum  number  of  customers  waiting  or 
being  served  equal  to  N,  the  total  number  of  UAVs  and  s  servers  (maintenance  paths. 
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each  with  service  rate  p).  Let  <p(x)  be  the  long  run  proportion  of  time  there  are  x  UAVs 
either  waiting  for  or  being  repaired. 

where  0(0)  is  determined  so  that  the  sum  is  one  and  N  is  the  total  number  of  UAVs. 

Approximate  the  long  run  fraction  of  time  during  which  there  is  at  least  one  UAV  on 
station  as 

«j(l)  =  (1-0(W))  a  <—  From  (5.7) 

long-run  fraction  of 
time  at  least  one  UAV 
is  up  from  queuing  model 


«i(2)=  (l-[«N)+*(iV-l)J  a 

v  V  '  / 

long-run  fraction  of  time  at  least  2  UAVs 
arc  not  waiting  for  or  being  repaired 

and  a  is  the  long  run  fraction  of  time  there  is  a  UAV  on  station  for  the  schedule  model  of 
Section  5. 

Numerical  Example 

Take  p  =£^+  x/(x+2T+^))p'  SuPPose  there  ^ 4  UAVs,  one  maintenance 
path  and  variable  ingress/egress  time.  For  the  parameter  values  of  the  previous  examples 

Fraction  of  Time  at  Least  One  UAV  is  On  Station 


Simulation  1 

.9^5% 

Confidence  . 
^Interval  ?  • 

:  v:-  Schedule  ^ 

! 

always  If 
•>.f  available  1| 

Schedule  Model 

at  least  one  UAV  i 
is  available 

;  Schedule  Mpdlel  ; 
time  Probability 

:a|least;t^l)i^Sj 

111  are  availahie  ::  II 

0.95-0.96 

0.96 

0.95 

0.91 

0.90-0.92 

0.91 

0.90 

0.86 

0.83-0.86 

0.87 

0.85 

0.80 

0.77-0.79 

0.82 

0.80 

0.75 

0.65-0.68 

0.76 

0.73 

0.64 

0.52-0.55 

0.71 

0.63 

0.52 
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APPENDIX  I 


B  is  the  time  between  takeoff  and  return  to  base.  Then 
f 2X  if X<T  (fails  on  flight  out) 

B  =  -jX+T  if  T  <X  <T  +  S  (fails  during  surveillance  time) 


(All) 


Thus, 


[2T  +  S  if  X>T  +  S  (fails  on  flight  back  or  does  not  fail) 

T 

E[B]  =  j2sAe~Asds 
0 

s 

+e-xrJ(s+2T)Ae~*sds 

0 

+e~X(r+S)(2  T+S) 

Let  A  be  the  time  on  station 

0  if  X  <  T  (Platform  fails  on  ingress) 

X-T  if  T <X  <T +S  (Platform fails  during  surveillance) 
S  ifX>T  +  S 


(AI.2) 


A- 


(AI.3) 


Thus, 


E[A]  =  0[l  -  *>-*■]+  e-*]  sXe^ds + e~x(-T+S)S 


-IT  1 


-XT 

=  e 


The  expected  length  of  a  cycle  is 


[I-,-] 


(AI.4) 


E[C]  =  E[B]+D+(l-e~x(2T+S))  - 


(AI.5) 


prob  of  failure 
during  cycle  expected 
repair  time 
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The  long-run  proportion  of  time  the  platform  is  performing  surveillance: 


n- 


m 

m 


-XT 

1  r,  -isi 

e 

H 

H- 

I 

■?[ 

l-e'w]+D  +  fl-e“(2r+5) 

J  A1 

J  L 

(AI.6) 
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APPENDIX  II 

A  Closed  Queueing  Network  Model 

Consider  a  closed  queuing  network  with  4  “service  centers”  or  states  labeled  as 
follows: 


Center 

Number  of  Servers 

D 

Down 

s  (number  of  maintenance  paths) 

I 

Ingress 

infinite 

S 

On  Station 

infinite 

E 

Egress 

infinite 

Consider  a  closed  migration  process  in  the  sense  of  Kelly  (1979).  The  state  of  the 
process  is  {(XD(0,  X/(f),  Xs(t),  XE(t));  t>  0}  where  XD(t)  is  the  number  of  platforms  down 
at  time  t;  X/(t)  is  the  number  traveling  to  the  region  at  time  t;  Xs(t )  is  the  number  on 
station  at  time  t,  and  XE(t)  is  the  number  traveling  from  the  region  at  time  t.  Assume 
{(Xo(f),  Xr(t),  Xs(t),  XE(t));  t  >  0};  is  a  continuous  time  Markov  chain  with 
P{ (Xj (t  +  h)  =  x j  - 1, Xj (t  +  h)  =  xj  +l,Xk(t  +  h)  =  xkfoTk*  i,  j\ Xe(t)  =  xt  for  all £} 

~  Xj \j(p  j  (Xj  +  o(h) 
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Put 


*di  =£[£>] +(l-e-A(2r+s))i 

A  Dj  =0  for  j  =  S,  E 

(pD{xD)  =  min(xD>s) 

Xjd  =  A 

Xis  =j 

A/y  =  0  for  j  =  E 

<Plixl)  =  xl 
A  5E  =B{A] 

A  $j  =  0  for  j  =  D,  I 

<Psixs)=xs 

Xed  =~ 

A  £j  =0  for  j  =  /,  S 
?£  (*£)  =  *£ 

Theorem  2.3  of  Kelly  (1979)  states  that 

lim  P  |(X£>  (/),  X/  (t),X5(t),X£(r))  = 

^n-rr2— 

r=l 

where  X  is  a  normalizing  constant,  0  =  (£>>  /,  S,  E }  and  §  is  the  unique  solution  to  the 
system  of  equations 
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=  £«**#  for  j  =  D,  I,  S,  E 
k  k 


Z«yi 


Numerical  Example 

Suppose  there  are  4  UAVs;  MTBMAF=25;  there  is  one  maintenance  path;  and  the 

other  parameters  are  as  before  with  the  ingress/egress  time  T  variable 

Long-run  proportion  of  time  at  least  1 UAV  is  on  station 

1  server 


T 

Limited 

Service 

Model 

Queuing 

Network 

Model 

Simulation 

1 

0.76 

0.79 

0.93 

2 

0.70 

0.74 

0.88 

3 

0.64 

0.68 

0.83 

4 

0.58 

0.61 

0.75 

5 

0.50 

0.54 

0.60 

6 

6.42 

0.45 

0.48 

If  there  are  4  servers 

Long-run  proportion  of  time  at  least  1  UAV  is  on  station 

4  servers 


1 

Independent 
Alternating 
Renewal  Process 
Model 

Limited 

Service 

Model 

Queuing 

Network 

Model 

Simulation 

1 

0.95 

0.96 

0.96 

0.95 

2 

0.91 

0.92 

0.93 

0.91 

3 

0.86 

0.87 

0.88 

0.86 

4 

0.79 

0.80 

0.82 

0.83 

5 

0.71 

0.71 

0.74 

0.78 

6 

0.61 

0.61 

0.64 

0.63 

In  the  present  case  the  multi-state  queuing  model  approach  appears  from  the  limited 
tabulation  to  be  biased  low  for  a  limited  number  of  repair  servers  (1).  This  is  not 
conclusive,  but  in  view  of  the  very  good  results  obtained  by  the  simpler  Markovian 


15 


model  of  Section  4  it  does  not  now  seem  desirable  to  use  the  present  approach.  Some 
further  investigation  may  be  worthwhile. 
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APPENDIX  m 

Display  of  Spreadsheet  Implementation  of  Analytic  Models 


ANALYTICAL  MODELS  FOR  MOBILE  SENSOR  (UAV)  COVERAGE  OF  A  REGION 

PARAMETER  VALUES 

1/Iam:mean  time  between  mission  affecting  failure: 

MAF(hours) 

25 

T :  ingress-egress  time 

3 

E:  total  endurance  time 

20 

alohau:  rate  of  unscheduled  maintenance  activities  (number 
per  hour) 

1/5 

1 /alphas:  time  between  scheduled  maintenance  activities 
(hours) 

50 

1/betau:  mean  time  to  repair  each  unscheduled  maintenance 
(hours) 

1.9 

1 /betas:  mean  time  to  repair  each  scheduled  maintenance 
(hours) 

7 

1/mu:  mean  time  to  repair  mission  affecting  failure 

2.4 

s:  number  of  maintenance  paths 

1 

N:  number  of  UAVs  (max  =*10) 

4 

S:  time  on  station  if  no  failures 

14 

1/betaa:administrative  and  logistic  down  time  (ALDT) 

0.5 

MOEs 

Long  run  proportion  of  time  at  least  1  UAV  is  on  station:  limited 
service  model 

0.617522645 

Long  run  proportion  of  time  1  UAV  is  on  station:  1  UAV  and 
alternating  renewal  process  model 

0.378147644 

For  Infinite  repair  capacity,  the  fraction  of  time  at  least  one 

UAV  is  on  station 

0.850462838 

CALCULATED  PARAMETERS 
exp  model  prob  of  failure  during  flight 

0.444444444 

lam:  MAF  failure  rate 

0.04 

prob  of  MAF  during  flight 

0.550671036 

prob  of  MAF  during  ingress 

0.113079563 

prob  of  MAF  during  surveillance 

0.428790936 

expected  time  of  flight 

15.16156427 

expected  repair  time  (Alt  Renewal) 

9.980959426  9.672902311  (limited  service 

expected  time  flying  UAV  is  on  station 

9.507586109 

proportion  of  flying  time  spent  on  station 

0.627084774 

gam:  rate  of  return  for  repair 

0.065956255 

X:  number  of  UAVs  down  long  run  prop  of 
time  number  of 

UAV  are  down 

0 

0.053634346 

1 

0.136872366 

2 

0.261968949 

3 

0.334266308 

4 

0.213258031 

5 

0 

6 

0 

7 

0 

8 

0 

9 

0 

10 

0 
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calcl 

1 

2.551953647 

1.913965235 

1.275976824 

0.637988412 

0 

-0.637988412 

-1.275976824 

-1.913965235 

-2.551953647 

-3.189942059 


calc2 

1 

2.551953647 

1.913965235 

1.275976824 

0.637988412 

0 

0 

0 

0 

0 

0 


ca!c3 

1 

2.551953647 

4.884350563 

6.232318117 

3.976146737 

0 

0 

0 

0 

0 

0 


ca!c4 

18.64476906 


calc5 

4 

3 

2 

1 

0 

0 

0 

0 

0 

0 

0 


calc6 

0.980660713 

0.948140259 

0.860934234 

0.627084774 

0 

0 

0 

0 

0 

0 

0 


calc7 

0.052597096 

0.1297742 

0.225538037 

0.209613312 

0 

0 

0 

0 

0 

0 

0 


E[repairtime]: 
Model  alt.  renewal: 


£[maintenance  time]+^1-e_^j— 


Model  limited  service: 


£[maintenance  time]  + 


A  1 

A-h-j-  /i 
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